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Методы распознавания образов и выделение 
кластеров при обработке файлов временной 
динамики температуры расплава 
кристаллизирующегося металла 


Рассматривается информационная модель теплофизических процессов при кристаллизации металла 
для дифференциального термоанализа (ДТА). Анализируется физический смысл высокочастотной 
составляющей ДТА-сигнала. Построен алгоритм отыскания кластеров в многомерном пространстве 
параметров эксперимента для последующего выявления структурных фаз в получаемом металлическом 
слитке. 


Введение 


Для различных задач металлургии, технологий литья металлов актуален автомати- 
зированный экспресс-анализ временной динамики температуры кристаллизирующегося 
расплава металла (динамический термоанализ — ДТА) [1]. Непосредственно получаемым 
файлом на первой стадии ДТА является цифровая последовательность показаний 
температуры термопары в пробнице с металлическим расплавом для совокупности 
равномерных отсчетов времени. 

Далее полученный файл отсчетов температуры обычно подвергается той или иной 
специальной процедуре сглаживания и/или фильтрации для получения некоторой 
достаточно гладкой кривой зависимости температуры от времени Т(®, которая затем 
подвергается однократному дифференцированию по времени и нормированию. 

Для получаемой при этом кривой (АТ(/а®/Т(@® = атТ(@)/а рассчитываются ее 
экстремальные точки и другие особенности, ассоциированные с началом или концом 
формирования основных фаз кристаллизирующего расплава в соответствии с файлом 
библиотеки файлов фаз многокомпонентных сплавов (БФМО) [2-4]. Как известно, данные 
о фазовом составе получаемого металла являются весьма востребованным компонентом 
информационной модели как самого процесса кристаллизации, так и всей совокупности 
предшествующих ему этапов подготовки и обработки анализируемого металлического 
расплава, а также степени адекватности итоговых данных ДТА [5] реальным процессам 
кристаллизации. 

Далее с целью изучения вклада явлений конвекционного переноса тепла в расплаве 
ставится задача на базе эксперимента выяснить меру зависимости показаний термопары в 
процессе остывания от места ее расположения внутри ограниченного объема пробницы с 
расплавом. Такой эксперимент направлен на выявление корректности учета неравномер- 
ности в пространстве и времени поля температур в кристаллизирующемся расплаве, 
находящемся в пробнице. Иными словами, анализируем степень однородности в 
пространстве процесса остывания и последующей кристаллизации расплава. 
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Получение температурных кривых 


Для задачи динамического термоанализа с участием автора [5] был проведен 
эксперимент на специальном моделирующем стенде «Мульти ДТА» по одновремен- 
ному снятию температурных кривых с 4 «почти одинаковых термопар» типа ХА, 
введенных через крышку сосуда в тестируемую жидкость симметрично (сосуд с водой, 
нагретой до точки кипения воды при давлении 1 Атм), как показано на рис. 1. 


Рисунок 1 — Схема расположения 4 термопар в сосуде с тестируемой жидкостью 


При этом получаемые данные 4 каналов в реальном времени поступают в 
компьютер через устройство ОЗВ М 9211 фирмы Манопа! плате. 

Предполагалось с целью извлечения необходимой информации на фоне шумов 
полученные при этом синхронные измерения четырех кривых охлаждения сопоставить 
между собой, используя результаты вычисления 6 парных корреляционных функций для 
4 температурных кривых, то есть для каждой пары каналов. В определенные моменты 
времени дважды (для { = 750 си{ = 1700 с) включался на 20 секунд дополнительный ис- 
точник тепла. 

Результаты прямых 4-канальных измерений поля температур приведены ниже на 
рис. 2. На рис. 3 приведены в качестве примера полученные в нашей работе [5] попарные 
корреляции для пар 1-2;2-—3;1-3;2—4. 


АА, о ра. 
Вр, ЧТО о х. 


0 200 400 500 800 1000 1 200 1 400 1 600 1 800 2000 2200 2400 2500 


Рисунок 2 — Результаты прямых 4-канальных измерений 
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Рисунок 3 — Некоторые попарные корреляции 
а), 6), в), г) для четырех кривых, показанных на рис. 2 


Результаты вычисления парных корреляционных функций приведены на рис. 3. 
Проанализировав рис. 3, можно сделать следующие выводы: 


1. «Шумовая» составляющая сигнала значительна: меньше 5 % при [Е (0,670]с и 


5...12 %1е (671,2400]с. 


2. В комплексе сигнал + шум доля шума скачком возрастает на участке { > 671 с. 

3. Из анализа парных корреляций сигналов 4-х каналов, показанных на рис. 2, 
следует, что подавляющую долю эксперимента (>95 %) сигналы в разных каналах 
коррелированны. Например для | и 2 каналов корреляционная функция превышает 0,5 
для 75 % времени и лежит в течение 10% времени эксперимента в области отрица- 
тельных значений корреляции (-0,2...—0,8). Подобная закономерность коррелирован- 
ности оказалась свойственна и парам 1 -— 3; 2 — 4; 3 - 4. Для пар 2 —-Зи 1 - 4 значения 
корреляционной функции не попадают в область отрицательных значений и лежат в 
пределах 0,6...1,0 в течение 62 % всего времени эксперимента. 

4. Из выводов | — 3 следует, что для понимания вклада реальных процессов 
конвективного теплообмена в охлаждающемся расплаве использование единственной 
термопары для большинства практических требований контроля качества получаемого 
слитка металла совершенно недостаточно. Необходимо проводить измерения парал- 
лельно значительным количеством однотипных термопар. 

5. При наличии данных М№хМ штук термопар получить данные измерения темпе- 
ратуры кристаллизирующегося расплава по МхМ каналам можно с помощью алгоритма, 
описанного в работе [6], рассчитать поле температур в приграничной (к внешней поверх- 
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ности) зоне изложницы. Точность такого метода расчета оценочно может быть получена 
методами работ [5], [6]. 

Это предположение означает, что высокочастотная составляющая ДТА-кривой (от 
нее мы избавляемся с помощью сглаживания сплайнами) несет важную информацию о 
ходе процесса кристаллизации конкретного образца сплава, которая в высокой степени 
индивидуальна для каждого конкретного образца кристаллизирующегося сплава. Это 
последнее предположение об информативности высокочастотной составляющей ДТА- 
кривой подтверждается результатами нашей работы [6] о результатах проведения 
эксперимента на специальном стенде «Мульти ДТА» для нескольких термопар 
одновременно как аналоговой модели теплофизических процессов в пробнице с 
кристаллизирующимся металлом. 

При исследовании (теоретическом и экспериментальном) движения турбулентных 
потоков жидкости или газа для обнаружения закономерностей измеряемых стохасти- 
ческих процессов в условиях быстрых пульсаций измеряемых параметров активно 
используют аппарат статистического анализа. Такой аппарат позволяет выявить 
некоторые важные закономерности для термобарических и кинематических параметров 
среды, усредненных по малым интервалам пространства и времени. Эти статистические 
методы включают, в частности, методы вычисления корреляционных и автокорреляцион- 
ных функций измеряемых параметров (поле скоростей, температур, давления и др.). 
В этом смысле, обнаруженные особенности поведения парных корреляционных функций, 
проявляющиеся на рис. 6, позволяют сделать предварительный вывод о возможности 
проявления в подобных экспериментах упорядоченных локализованных 3З-мерных 
вихревых структур, которые относительно медленно распространяются в тестируемом 
объеме остывающей жидкости, перенося с собой запасенную в данном локальном объеме 
соответствующую порцию кинетической энергии, которая в конечном итоге со временем 
диссипируется во всем объеме жидкости. Такие объекты относятся к категории так 
называемых ползущих течений и в нашей работе [6] условно названы «блуждающими 
вихрями». 

Наличие системы «блуждающих вихрей» может давать значительный вклад 
в снижение градиента температур в охлаждающейся жидкости. Это может быть исполь- 
зовано как дополнительный механизм управления процессами кристаллизации расплавов 
металлов, полимеров и других материалов, получаемых в металлургии, органической и 
неорганической химии и физической химии. Такой вклад альтернативен вкладу в 
тепломассоперенос обычной теплопроводности, которая реализуется за счет механизма 
чисто молекулярной тепловой диффузии. 


Очистка базы данных температурных кривых 


Для построения системы автоматического определения характеристик алюминие- 
вых сплавов методом дифференциального термического анализа необходимо иметь базу 
данных температурных кривых определенного качества, о котором говорится далее. При 
записи температурной кривой с помощью пробницы с термопарой особые требования 
предъявляются к спаю термопары. Накопленная ранее база данных температурных 
кривых не может быть полностью использована в качестве обучающей выборки, так как 
содержит термограммы, по которым не проводились дублирующие методы анализа. 
Формирование обучающей выборки является основной целью данного исследования. 

Процесс кристаллизации алюминия в пробнице протекает за время, определяемое, 
главным образом, фазовым составом. Поэтому в обучающую выборку выбираются 
кривые определенной длины. Длина температурной кривой зависит от температуры 
начала записи (жидкий сплав 700°С) и температуры конца записи (твердый сплав 400°С). 
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Используемый алгоритм ДТА предполагает вычисление кривой темпа охлаждения. 
Исходя из этого, возникает ограничение на гладкость исходной кривой. В работе 
проводится сглаживание кривых температуры и темпа охлаждения с помощью сплайнов. 
В большинстве случаев значения темпа охлаждения попадают в определенный диапазон. 
В связи с этим отбраковываются кривые с большой разницей максимального и 
минимального значений темпа охлаждения. 
На следующих рис. 4 -— 6 показано, что подавляющая часть кривых темпа 
охлаждения имеют небольшие значения минимального тш_т и максимального тах_т 


темпов охлаждения. На рисунках ось абсцисс соответствует минимальному значению 
темпа охлаждения, а ось ординат — максимальному. 


Е РА И СЕ ЕЕ 


ен 


300$ ---------- < к------------- =. ь----- 
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А 200 1 000 -800 -600 -400 -200 0 


Рисунок 4 — Зависимость тах_ т(т_ т) 
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Рисунок 5 — Зависимость тах_ (та _ г) (увеличено в 10 раз) 


-40 -35 -30 -25 -20 45 0 -5 0 
Рисунок 6 — Зависимость тах_ (та _ г) в окрестности основного кластера 


На рис. 7, 8 показано, что обучающую выборку можно составить из 77 % или 85 % 
температурных кривых, для которых разница максимального и минимального 
темпов охлаждения лежит в пределах 20 единиц. 
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Рисунок 8 — Доля удовлетворительных кривых после сглаживания 


При этом 85 % кривых составляют кластер с координатами (-7,486; 1,004) и 
дисперсией (3,082; 1,144). 


Выводы 


Сопоставление гистограммы до и после сглаживания приводит к предположению 
о том, что сглаживание сплайнами, как оказывается, уменьшает индивидуальность 
отдельных ДТА-кривых (формально за счет «мягкого игнорирования высокочастотного 
шума /, + АЁ Гц»). Это предположение, по сути, означает, что высокочастотная состав- 


ляющая ДТА-кривой, от которой мы избавляемся с помощью сглаживания сплайнами, 
несет важную информацию о ходе процесса кристаллизации конкретного образца 
сплава, которая в высокой степени индивидуальна для каждого конкретного образца 
кристаллизирующегося сплава. Последнее предположение об информативности высоко- 
частотной составляющей ДТА-кривой в значительной степени подтверждается резуль- 
татами нашей работы [6] о результатах проведения эксперимента на специальном стенде 
«Мульти ДТА» для нескольких термопар одновременно как аналоговой модели 
теплофизических процессов в пробнице с кристаллизирующимся металлом. 
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ТА. Варава 

Методи розшзнавання образйв 1 видлення кластер1в при обробщ файлв часово! динам ки 
температури розплаву металу, який кристалзуеться 

Розглядаеться 1нформащийна модель теплофзичних процесв при кристалзаци металу для диференщального 
термоаналзу (ДТА). Аналвуеться флзичний змст високочастотно! складово! ДТА-сигналу. Побудований 
алгоритм в1дшукування кластер1в у багатом1рному простор! параметрав експерименту для наступного 
виявлення структурних фаз в одержуваному металевому злитку. 
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